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On the Nature of the Diffraction Figures due to the Heliometer. 
By Philip Francis Everitt, University College, London. 

(Communicated by Prof. Karl Pearson, F.R.S. — Received September 11, — 

Read November 25, 1909.) 

[Plate 3.] 

1. In investigating the influence of personal equation in connection with 
double star measurements, the actual values of distance and position angle 
must be known, and hence the only possible method is to use an artificial 
star. The method adopted in the Students' Observatory at University 
College, in order to obtain a known separation and direction, is to use a 
weak heliometer lens in front of the equatorial and a spark between cadmium 
electrodes some distance away. The direction of the line joining the two 
" star "• images and the distance between them can be measured by a divided 
circle and micrometer screw respectively, while the intensity of the stars can 
be altered by using various stops. It was found that when a moderately 
small stop was used, the diffraction images became quite perceptible, and 
hence, before any further progress could be made in the determination of 
personal equation, it was necessary to investigate these diffraction images. 

2. On consulting the literature of the subject, I found that the earliest 
papers dealing with the diffraction figure of the heliometer are by Bessel. He 
observed (a) that two opposite brushes of light emerged from the central disc 
perpendicularly to the line of separation ; (b) that when the halves of the 
lens were separated, the disc was lengthened in the direction perpendicular 
to the line of separation ; and (c) that the brushes of light were of equal 
length when the halves were united ; but that, when the halves were 
separated, the brush of light on the same side of the image as the half -lens 
forming it was longer than the one on the opposite side ; and he attributed 
the lengthening of the disc formed by the half-lens to the fact that the lens 
was corrected for spherical aberration for the whole lens, and not for the 
two halves separately. 

This explanation does not, I think, require further examination. He also 
observed that there was an instability of superposition of the two images, 
and this was confirmed by Hansen and Gauss. In 1882 Bruns treated the 
diffraction figure of the semicircular aperture analytically, and gave formulae 
for the intensity of the light at any point. These formulae demand, how- 
ever, a very large amount of arithmetical work before they can be of practical 
use. Bruns stated further that the task of computing the distribution of 
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light in the image had been undertaken by Schnauder. The only other paper 
on the subject was by Schemer and Hirayama. It consisted chiefly of plates 
of photographs of diffraction images formed by various apertures, including 
the semicircle, and a brief summary of the paper by Bruns already 
mentioned. There are also notes and a definite statement that the calcu- 
lations initiated by Prof. Bruns had not been completed. 

3. I decided to attack anew the problem of computing the diffraction 
figure, hoping by the help of modern mechanical integrators to be more 
successful. The formulae developed by Bruns give the light intensity, H, in 
terms of two components and S, such that 

H = C 2 + S 2 , 
where 

C = 7T7 2 ^4 , S = 2v 2 - f— — — cosw )du, and ^ = Scos^: 

8 J Q v\ v I 

8 and e being the ray and angle polar co-ordinates of any point of the 
diffraction image referred to the centre of the image as pole, and a line 
parallel to the diameter of the aperture as initial line, and 7, a constant such 
that, if the intensity at the centre of the image be taken as unity, we must 

put 7T7 2 =2. 

The values of C are obtained without difficulty from any table of Bessel 
functions for a somewhat limited range of values of 8, but the values of S 



depend on the evaluation of the integral I 



' e 1 /sin v x 7 

v\ V 



Three methods suggest themselves, namely (1) to make this integral depend 
on other tabulated integrals, (2) to expand in series, and (3) to evaluate by 
mechanical integration. It was not found possible to use the first method, 
and so only methods (2) and (3) were available. Since v = 8 cos u, it is clear 
that the choice of method depends largely upon the range of values of 8 that 
are to be considered. JSTow 8 in the portion of the diffraction figure actually 
drawn ranges from to 16*47. For the central region from to 3*83 
expansion by Taylor's theorem was used, and for the remaining portion, where 
expansions of the integral converged with extreme slowness, mechanical 
integration was adopted. 

4. Evaluation by expansion, using Taylor's theorem. With the usual 
notation, let /(0) be the value of the function I and f 1 (0) the value of 
the nth derived function when 8 is put equal to after obtaining the derived 
function. 



Let 



e 1 /sin v 1 7 
cos v du, 

v \ V 



Y 2 
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then 









IS 2 = 



" e sin v 



w 



du 



COS V 



du 



Jo 



V 



sin (B cos it) du 



o 



S 2 cos* u 



cos (S cos u), du . 
S cos u 



sin (S cos w) du f e S cos (S cos u) du . 



o 



cos 2 u 



d_ 

dS 

dJ 



/t£2\ _ f 6 cos P (S cos u) du 
J COS u 



cosu 
cos*(S cos %) du 



COS 16 



-+• | S sin (S cos u) du, 
o 



or S — + 21 = sin (S cos iA cZ% ; 
dS Jo 



.'. /(0) = 0. 

Differentiating again, 

cos (S cos u) cos % ^ ; 



o dl . £ d 2 l _ 
Proceeding, 4 



o 



•• /'(0) = *J 



= -J I cos u du. 

o 



cos 2 u sin (S cos %) du y 







whence /"(0) = 0. 

As we proceed further we shall always, in the case of the derived functions 
of even order, have on the right-hand side, as factors under the integral sign, 
an even power of cos u and also a sin (S cos u), consequently f 2n (0) = 0. 

For the odd derived function, the law is best expressed as follows:- — 



(2n + 3)/ 2w+1 (S).+ Sf 2n+2 (S) = ( - If I cos 2w+1 ^ cos (S cos u) du, 

o 



whence 



(2^ + 3)/ 2 » +1 (0) = (-l)M cos 2 ^ +1 

Jo 



^6 du. 



Six terms in the expansion, i.e. proceeding to f xi , were found to be 
sufficient for the range of fairly small S's dealt with by this method. 

The values of the integrals cos 2n+1 udu were obtained by expressing 

Jo 

the powers of cos u in terms of cosines of multiples of u and then integrating. 
The values so obtained were — 
/i(0) =isin€, 

fm(Q\ - ! /sin 3 e 3 sin/ 

/v/n\ — i / sin 5 e 5 sin 3 e 5 sin e" 
7 W -T^-gQ™+ 48 + 8 



/ vii (0) = -i{ 



sin7e , 7 sin 5 e , 7 sin 3 e , 35 sine 



448 



4- 



320 






64 



64 
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1 fsin9 e 9 sin 7 e 9 sin 5 e 7 sin 3 e 63 sine 



/ix (0) = JL J i 

w l ] 1 2304 1792 320 64 128 



+ 



+ 



■ + 



}• 



/xifft\ = _JL 



(0) 



1 3 



sin 11 e . 11 sin 9 e 55 sin 7 e ,33 sin 5 e , 55 sin 3 e 

I — ~ ■■>., ^ + _ ZZZ nn — T T7T^~ i I 



11264 



9216 



7168 



1024 



512 
, 231 sin 



1 sin el 
512 J ' 



In using this method it was found most convenient to calculate the 

S n 1 

general term ~ f n (0) as the product of B n and — f n (0), so that on one sheet 

were calculated the values of the latter function for values of e by intervals 
of 5° from to J tt, giving coefficients which, when multiplied by the 
respective values of B n and summed, provided the values of I. 

The calculations were rendered less tedious by the use of a Brunsviga 
calculator. 

Table I.*— Calculation of — f n (0). 





5°. 


10°. 


15°. 




sin 6 


-0871557 


-1736482 


-2588190 




sin3g 


'2588190 


-5000000 


-7071068 




sin 5s 


'4226183 


-7660444 


-9659258 




sin 7s 


0-5735764 


-9396926 


-9659258 




sin 9 6 i 


-7071068 


1 -ooooooo 


-7071068 




sin lit 


-8191520 


-9396926 


-2588190 


tit sins 


-0290519 


-0578827 


-0862730 




sin 3s 


0*0215683 


-0416667 


-0589256 


12 




3 sin £ 
■— 4 — 


0*0653668 


-1302362 


-1941143 




sum 


-0869351 


-1719029 


'2530399 




-io sum 


-00289784 


-005730097 


-008434663 




sin oe 
~80~' 


0*00528273 


0-00957556 


-01207407 


■ 


5 sin 3e 

48 


-02696031 


-05208333 


-07365696 




5 sin f 


0-05447231 


-10853013 


-16176188 


8 




sum 


-08671535 


-17018902 


-24749291 




¥Io sum 


-0001032326 


-0002026060 


-00029463442 



* This table and those on pp. 307 and 308 were calculated to more figures than are requisite for 
the final results, as, a priori, the degree of approximation required could not be determined. 
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Table I. — continued. 



sin 7s 

~~448" 

7 sin5e 

7 sin 3e 



64 
35 sin e 



64 



5°. 



10°. 



15 c 



'00128030 
-00924478 
'02830833 
-04766327 



-00209753 
-01675722 
-05468750 
-09496386 



-00215608 
-02112963 
-07733981 
'14154164 



sum 

sum 

45360 



sin9e 

"2304 

9 sin 7e 

1792" 
9 sin 5e 

"~~320~~ 

7 sin 3e 

63 sin e 
__ _____ 



sum 

sum 



3991680 



sin lie 

11264" 
11 sin 9e 

~~9_l6~" 

55 sin 7e 

7168 
33 sin 5e 

"1024~ 

55 sin 3s 
__ ______ 

231 sine 
_____ 



sum 
sum 



518918400 



-08649668 

-01906893 x 10~ 4 



-16850611 

'03714861 x 10~ 4 



'24216716 

-05338782 x 10~ 4 



-00030690 
-00288069 
-01188614 
-02830833 
-04289695 



-00043403 
-00471944 
-02154500 
-05468750 
-08546747 



-00030690 
-00485119 
-02716666 
-07733981 
-12738748 



-08627901 

'02161471 x 10" 6 



-16685344 

-04180030 x 10 ~ 6 



-23705204 

-05938653 x 10~~ 6 



-00007272 
'00084399 
-00440105 
-01361954 
-02780282 
-03932220 



-00008342 
'00119358 
-00721025 
-02468698 
'05371094 
-07834518 



-00002298 
'00084399 
-00741054 
0*03112847 
-07595874 
-11677185 



-08606232 

-16584942 x 10~ 9 



-16523035 

-31841297 x lO" 9 



-23213657 

-44734696 x KT 9 
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Table II. — Calculation of I. 







5°. 


10°. 


15°. 


5 = 2 0— 

5 

8 3 


2 

8 

32 

128 

512 


0-0581038 
-0-0231827 

-0033034 
-0 -0002441 

-0000111 


-1157654 
-0-0458408 

-0064834 
-0-0004755 

-0000214 


-1725460 
-0-0674773 

-0094283 
-0-0006834 

-0000304 


5 5 

57 


$> 








-0379915 


-0759539 


0-1138440 


d = 2 -4— 
d 


2-4 
13 -824 
79 -6264 
458 -647142 
2641 -807538 
15216 -81137 


-0697246 
-0-0400597 

-0082200 
-0-0008746 

-0000571 
-0-0000025 


-1389185 
-0-0792129 

-0161328 
-0-0017038 

0001104 
-0-0000048 


-2070552 
-0-1166003 

-0234606 
-0-0024486 

-0001569 
- -0000068 


P 


8* 

*7 


p 


c 11 








-0370649 


-0742402 


'1116165 


d - 2 -7— 
£ 


2-7 

19 -683 
143 -48907 
1046-0353 
7625 -5973 
55590 -6043 


-078 1401 
-0-0570382 

-0148127 
-0-0019947 

-0001648 
-0-0000092 


-1562832 
-0 -1127855 

0-0290717 
-0-0038859 

-0003187 
-0-0000177 


-2329371 
-0*1660195 

-0422768 
-0-0055846 

-0004529 
-0-0000249 


P 

5 5 

£7 


£9 

d U 








0-0343755 


-0689845 


-1040378 



5. Evaluation by mechanical integration. For the values of S between 
3*83 and 16*47 mechanical integration was employed. Here it was found 
desirable to alter the form of the expression to be integrated, so that it 
would be less cumbersome in the work of computing the points on the 
curve before plotting. 



Putting 



1 /sin v \ 

z = cos V , 

V \ V J 



then I = 



z du 







and if we plot z against u, the area up to any value e of u will be the value 

of I for that e. 

Put v = S cos u = tan <£, 

17 sin 1? \ sin-y cos 6 — cosi? sin<f> 
; _ cosi? = j~ r.. 

v \tan <j> 1 v sin <£ 

therefore z = \ ""**? / , 

v sin (j> 

and this is a form eminently suited to logarithmic calculation. 



then 



z 
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The work was arranged in parallel columns, each column for one value of 
u and each sheet containing all the work for one value of S. For the 
preliminary calculations, values of u were taken at intervals of 5°, the 
values of z so obtained were plotted on squared paper, and the maxima and 
minima (which will be separately discussed) inserted. From this graph it 
was evident what further points were required in order to draw the curve 
of . z accurately. This rough plotting on squared paper served as a very 
useful check against errors in the work of calculation. The arrangement is 
shown in Table III. 

The positions of the maxima and minima were calculated by a different 
method. 

Table III. — Calculation of z. 
8 = 11-61986. 





0°. 


5°. 


10°. 


15°. 


Max. 


log 5 

log cos u 


1 -0652009 


1 -0652009 
9 -9983442 


1 -0652009 
9 -9933515 


1 -0652009 
9 -9849438 


1 -0652009 
9 -2531911 


log v = log tan <j) 


85" 4' 52" -6 


L -0635451 
85° 3' 45" -3 


1 -0585524 
85° 0' 20" -8 


1 -0501447 

84° 54' 31" -2 


-3183920 


, v in radians 
v in degrees 


11 -619860 
665° 46' 8" *2 


11 -575616 
663° 14' 2" -2 


11 -443330 
655° 39' 16" -1 


11 -223922 
643° 5' 0"'l 




v —<{> 


580° 41' 15" -6 


578° 10' 16" -9 


570° 38' 55" -3 


558° 10' 28" -9 




log sin <p 


9 -9983977 


9 '9983855 


9 -9983480 


9 -9982831 


o 


log ( 1 ) 

\sin<£>/ 


0-0016023 


0-0016115 


-0016520 


-0017169 


00 

zp 

o 
1.^ 


log sin (v — </>) # 


8 -9347991 


8 -9364549 


8 -9414476 


8 -9498553 


II 


9 -8142034^ 


9 -7909994^ 


9 -7073766^ 


9 -4940367^ 


?5 


log z 

z 


8 -7506048^ 
-0-056312 


8 -7290688^ 
-0-053588 


8 -6504762^ 
-0-044717 


8 -4456089w 
-0-027900 





* n in this table denotes that the numerical quantity of which the logarithm is given occurs 
with a negative sign. 



Since 



1 /sin v \ 1 <s 

% = _ cos v and v = o cos u, 

v \ v 



then 



dz dz dv dz , ^ . N 
— = — . — = _ . ( — bsmu). 

du dv du dv 



Hence, for any value of S, not zero, u not being zero, dzjdu and dzjdv have their 
zero values simultaneously ; consequently, instead of considering dz/du = 0, 
we may take dzjdv = 0. We have then — 





dz ■ 2 . , 2 , 1 . 

= — -- sin v + - cos v -f - sm v, 



dv 



v s 



w 



V 
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. • /2 1\ 2 eos^ . , 2v , 
or sm^ -— - = 5—, %.e. tan / y = -- 5 , 

\v 6 v J , w* Z — v^ 

and the roots of this equation give the values of v, and therefore of u, for 
which z is a maximum or minimum. 

The roots were found by successive approximation. 

Put v — Vo + e, 

then tan (v + e) = y y -^ . 

2 — (% + e) 2 

Now expand the left-hand side by Taylor's theorem and the right by the 
binomial, neglecting second and higher powers of e. 

Then 

4. , e 2lv + e) 2v , f 2 , 4 V 



cos 2 v 2 — ^o 2 — 2 ^ 2 — v 2 12— ^ 2 (2 — v 2 ) 2 



Therefore , 2 v 



2— ^o 2 2 — ^o 2 f, 2v 



e = ■ ■ — _ = -^ -< tan % 



< tan 



4^o 2 _ o3rt . ^o 2 L 2-w 



+ 77; — ^Tno — see^o 



2-V (2-^ 2 ) 2 



= J(l_ 2^ C ot, 
v \ 2—v 2 



This formula was not very often required, as is easily seen from the 
following considerations. Counting the maxima and minima from u = \ir 
backwards to u = 0, there is only a limited number of them for any given 
value of 8 ; moreover, these maxima and minima occur with precisely the 
;same values of v and z, when 8 takes up a new value, the only change with 
regard to them being that u changes its value in accordance with v = 8 cos u. 
The values of the maxima and minima may be calculated by the same method 
as the other values of z, or more shortly as follows : — 

2v 



tan# = 



2-v 2 
2 



v. 



therefore % = -.(— — ^ 

All the data required for the plotting of the curves of z for the different 
values of 8 being now obtained, the curves were plotted on large, thin, 
smooth-surfaced* cards in rectangular co-ordinates by means of a Coradi 
-co-ordinatograph, and the curves drawn in with the weights and splines 
used in naval architecture. The integration was then performed by 
means of a large Coradi integraph, which draws the graph of the integrated 
function. The cards were again ' placed under the co-ordinatograph, and 
* Because both co-ordinatograph and integraph work best on such cards. 
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the ordinates ruled for successive values of e, increasing by steps of 2°. 
The photograph of such a card is shown on a much reduced scale in 
Diagram 1. 



ED*ttttg'.' 



PD-4868f 





Ordinates of Pi h B' from sloping- base lines x P. D. = inCegral 
to corresponding- points of primitive curves A h B respectively. 



Diagram 1. — Mechanical Integration of the I-Integral. 
(Reduced to 1/6 linear of original drawing.) 

The values of 8 first taken were those giving the maxima and minima of 
the ring system of the circular aperture, and afterwards a sufficient number 
of intermediate values were taken to give an accurate curve when H was 
plotted against 8. The values of the ordinates of the integrated curve were 
scaled off and tabulated. For the angles 10°, 20° . . . 80°, 90° and 8 — 3*83 . . . 
the values of I were obtained by both methods, and these values were 
compared. The values are given in Table IV. 

Table IV. 





10°. 


20°. 


30°. 


40°. 


50°. | 60°. 


70°. 


80°. 90°. 


Direct calcula- 
tion 

By mechanical 
integration 


-0283 
'0283 


-0617 
-0616 


-1046 
'1046 


-1603 
-1603 


-2285 
'2290 


'3039 
-3044 


-3749 
'3752 


'4266 ; '4455 

i 

-4272 ! *4460 



These show a maximum difference of 1 in 457, or about 2 per mille. The 
greater the value of 8, however, the greater the inaccuracy of the mechanical 
integration, but the increase is slow ; but, on the other hand, we have also a. 
less accurate result from the less convergent series. 
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6. The evaluation and representation of H. 

The value of H was obtained by combining the two parts C and S. Putting 

2 4 

7T7 2 = 2, we have C = J\ (S) x -k, and S = — I. This part of the work was also 

carried out in tabular form, one sheet being used for each value of S. The 
values of H had now been obtained for a meshwork of points, but the bare 
figures did not help one to realise the appearance of the diffraction image. 
Eepresentation by means of a model was unfortunately out of the question^ 
as the ratio of the intensities in different regions of the figure is too great. 
It was found far more useful to construct a diagram showing contour lines 
or lines of equal light intensity. In order to do this, interpolation curves 
had to be drawn for every value of e, plotting H against 8. It was found 
necessary to draw, each curve in several overlapping sections, using different 

Major axis 



Intensity at centre- 1. 




000^ Minor axis and 
base of semicircle 



Diagram 2. — Contour Lines of Diffraction Figure. 

1 
(Keduced to — linear of original drawing.) 
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scales for each, in consequence of the large range in value of H. These curves 

were drawn upon squared paper, which for this purpose was found sufficiently 

accurate. The values of 8 } at which H took up any desired value, were then 

tabulated for the intervals of e adopted, and from this table the contour 

diagram was constructed. This Diagram (2) was plotted on a large scale, and 

then photographically reduced. The central intensity being unity, the values 

of H used were 0*9,, 0-5,. 0% 0*05, 0*02, 0*015, 0*006, 0*004, 0*002, 0*001, 

0*0005, and 0*0002. The contour diagram shows all the features seen when 

the diffraction image is examined visually, as is evident on comparison with 

the photograph taken for that purpose. 

7. Having obtained the diffraction figure, up to a distance of 8 = 16*47, 

two questions at once presented themselves : (a) what is the accuracy of the 

result, and (&) is it desirable to extend our knowledge any further. Dealing 

with the latter question first, it certainly does seem desirable to trace the 

-fluctuations of light along the bright ray which forms the major axis of the 

figure ; and this, being carried out by direct calculation, will serve as a test 

of the accuracy of the mechanical integration at the part of the range where 

the accuracy so obtainable is smallest. The formulae hitherto used will not 

serve when 8 is, as now, large. Brims gives as an alternative expression 

for S, 

2rv 2 r 4 T /CA . . 8 T ,*v . . ,12 



osma 



- — : J 2 (8) sin 2 a + 7r~p ^4 (*5) sin 4 a -j- zA^ Je (S) sin 6 « + . . . 
1.3 w 3.5 w 0.7 



where . & = c— |-7r. 

Now, for the major axis of the figure a = 0. 

Hence we have 

"4.2 T ,* x . 8.4 T ,*, , 12.6 



8 



J2(3) + ^|J 4 (8) + J ^fJ6(S) + ... 



1.3 v ' 3.5 v 7 5.7 
Putting in the value of y 2 , and taking a factor 8 out. of the bracket, we have 

q _ 32 g n 2 J 2n (5) 



which may also be written 



7r§7 4w 2 — 1 



8 f"£ -r ^x . 3 3^8 



7r8 l 2 4n 2 — 1 

The values of H were then calculated for 8 = 10, 14, 15, . . . 24 and 
plotted; intermediate values of H required for splining were obtained by 
interpolation, using a parabola of the third order, this being far quicker than 
interpolating for each J 2n (&) required. This method served up to 8 = 24, 
at which argument the tables of J's at present end. The values of H for 
€ = \ic were then tabulated and the graph of H against S (e = \nr) drawn. 
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The value of H was obtained by both methods in four cases and the results 
are given in Table V. jit. is at once evident that the mechanical integration 
gives accurate results, the maximum error barely exceeding one unit in the 
fourth place of decimals. 

Table V. 



5. 


io-o. 


14 -0. j 15 -0. 


16-0. 


Mechanical integration '] 
Direct calculation 


-0318 
-03189 


'0049 
-00493 


-0085 
0-00850 


-0104 
'01056 



The above method, as I have said, was available up to S = 24, beyond 

which the J's are not yet tabulated. As this range does not carry us beyond 

the third bright spot, and it was desired to go further than this, another method 

was necessary. After several futile attempts, such a method was found and 

has proved itself very useful. It would, of course, have been possible to have 

computed all the J's required (a very large number), but the labour would 

have been excessive and the following method was adopted to lessen the 

work: — 

sin (x sin <£) = 2 { J\ (x) sin <£ + J 3 (x) sin 3 <£+... }, 

I sin (x sin <£) d<j> = — 2 {Ji (x) cos <£ + •£■ J3 (#) cos 3 <£ + . . . }, 



f7T 



sin (x sin <£) d<j> = 2 { Ji (x) + 1 J 3 (x) + ^ J 5 (x) + . . . }. 
Differentiate both sides with regard to x — 

cos (x sin <£) sin cf> d(j> = 2 { JY (x) + -i- J 3 ' (x) -f £. J 5 ' (x) -f . . . } ? 



C( 

Jo 



but 
therefore 



J n '{x) = 1 {J n -l(%) — Jn+l(%)}, 



f 



*ir 



cos (x sin <£) d (cos <£) = { —Jo (#) — 1 J2 (#)— -5- J4 (#)—■ 7- Je (#) — • • • 



+ J2(x) + $Jt(x)+-$J 6 (x)+...} 

/ 1 1 



+ 



\2ra-l 2^+1, 



J&i («)+... 



00 



= -J (^)+2 



2 J 2?i (#) 



Therefore 



r 4^ 2 -l * 
^4^a-Ti='iJo(«) + i 



cos (# sin <£)*rf (cos <£).• 
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Substituting this value in 



we have 

or, putting z = cos <£, 



S = 



s = 



8 



7rS 

_8_ 

77*8 



Jo (8), g J2tt(S) l 



2 




i4:n 2 -lj ' 
cos (8 sin (/>) rf (cos <£) 



S = 



7T< 




ri 



o 



cos (8^/(1—$?) dz) 



The determination of S resolves itself, therefore, into the evaluation of 
this integral. The method adopted may be called a semi-graphic one and is 
identical with that used independently by J. H. Shaxby and described by 
him in a paper read before the Eoyal Society on April 22, 1909. 

The method had been found, tested and definitely adopted for this work at 
the end of the Lent term, and before I was aware by Mr. Shaxby ? s paper 
that he also had discovered its advantages. The method depends on the 
property possessed by this class of curves that, after cutting the axis a third 
time, the areas of all the other sections are almost identical with the areas of 
sine curves of the same altitude on the bases formed by the various points 
of intersection of the original curve with the axis. This is seen from the 
values given for 8 = 33 and S = 40 in Table VI. 

Table VI. — Areas of Sections. 





5 = 


33. 


5 = 


40. 


Calculated. 


Planimeter. 


Calculated. 


Planimeter. 


1st partial section 


i 

-2538 
-1024 


-oooi 

'^llo 
-1016 


'1522 

"0868 


'0873 
-1469 
-0862 




22x10. , , j j ......... 


3rd complete section. 


0-07L7 


-0717 


-0638 


'0636 


4th „ „ 


'0543 


'0545 


'0501 


-0501 


5th „ „ 


-0423 


-0424 


-0405 


-0404 




-0329 


-0329 


-0329 


-0331 


7th „ „ 


'0250 


-0249 


-0268 


-0269 


8th „ „ 


'0181 


-0181 


-0214 


-0216 


9th „ „ 


6 '0118 


-0117 


'0165 


-0164 


10th „ „ 


-0058 


*0059 


-0121 


'0120 


XJLuH. ,} ,, ......... 


-0007 


*0007 


'0080 


0-0080 


12th J5 „ 


• — • 


— , 


-0039 


'0036 


13th „ , 


! 


~ 


-0005 


*0005 



The horizontal line drawn through the above table marks where the curve 
cuts the axis for the third time. Beyond this point the calculated areas are 
practically identical with the planimetered areas of the curve. 
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The curve of y = cos B^/(l— z 2 ) was therefore traced from z = up to 
the point where it cut the axis for the third time and the area measured by 
a Coradi's Compensating Planimeter ; the other points of intersection were 
calculated, the areas of sine curves on the bases so obtained computed and 
the whole area of the curve or value of the integral obtained by addition. 

In dealing with graphs of this curve y = cos S(l — z 2 )%, the usual practice 
of calculating ordinates for arbitrarily assigned abscissas was reversed ; the 
ordinates were arbitrarily fixed such that they would give equal increments 
of abscissa if the curve were a sine curve and the corresponding abscissae 
were calculated for the actual curve ; this method giving the points of inter- 
section directly and also dividing the calculated points fairly evenly along 
the curve to be traced. 

Accordingly a sufficient number of values of x were taken and the 
ordinates were cos x. 

Then z the abscissa of our actual curve was given by 

'2 7rn-{-x\ 2 



z 



(M 



}' 



n being an integer and ranging from upwards so long as z remains real. 

Table VII. 





a. 


b. 


c. 


1 — --....I... ... .... . ■■ -.1— ..^ 


2irn 

X 


31 '41593 
1 -57080 


31 -41593 
1 -S9626 


31 -41593 
1 -22173 


31 -41593 ; 
1 -04720 

i 


2irn + x 


32 -98673 


32 -81219 


32 -63766 


i 

! 

32 -46313 


log (2irn + x) 

Ml) 


1 -51834 
8 -39704 


1 -51604 
8 -39794 


1 -51372 
8 -39794 


1 -51139 
8 -39794 


21qg( 2, "; +a, )-logJ 


9 -91628 
9-83256 


9-91398 
9 -82796 


9-91166 
9 -82332 


9 "90933 
9 -81866 


B = log (1) -log 5 

A 


0-16744 
9*67247 


-17204 
9 -68670 


-17668 
9 -70074 


-18134 
9 71453 


log z 2 = log b + A 
logs 


9 -50503 
9 -75251 


9 -51466 
9 -75733 


9 -52406 
9-76203 


9 -53319 
9 -76660 


z 


-56560 


-57191 


-57814 


-58425 
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In the numerical work of computing z (shown in Table VII) use was made 
of the tables of " Additions- und Subtraktions-Logarithmen," given in 
Bremiker's five-figure logarithmic tables, which are specially suited to the 
type of formula in question. The specimen in Table VII is written out at 
full length, but in the actual work the starting point was of course 
log (2 it n -f x) and the line B was omitted. 

The small letters over the columns indicate to which of the arbitrarily 
fixed ordinates the abscissa £ corresponds. 

The curve was, as before, plotted on cards with the co-ordinatograph and 
drawn in by the help of a spline. The complete curve for S = 40 is 
reproduced in^Diagram 3. 




Diagram 3. — Graph of y = cos 8(1 — z 2 )i for d = 40, to show approach of areas to those of 

sine curve. (Beduced to 1/4 linear of original drawing.) 

The planimetered areas J when added to the residual areas calculated by 
fitting halves of sine curves, gave the required integral, and consequently S. 
The value of was found from the values of Ji (8) as before; but these 
could not now be taken from tables of Bessel functions, and therefore were 
computed from the well-known semi-convergent series.* The accuracy of 
the values of H so obtained is very good indeed, as is shown by comparing 



* Gray and Mathews, 'Treatise on Bessel Functions,' p. 40. 
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the value obtained by the above method and that by summation of J's for 
the ease of S = 24, where the two methods link up — 

By summation of J's H = 0*003478 

By semi-graphic method ..................... H = 0*003479 

This unusual accuracy is due to the fact that the integral evaluated is small 
and is subtracted from unity, thus decreasing the percentage error, which 
here is only 1/30 per cent. 

The values of H have been tabulated for the cases of the two principal 
axes of the diffraction figure, e = and e = J 77*, and from these have been 
drawn the curves of intensity of illumination along these axes (Diagram 4). 
Of course the intensity along the axis e = is identical with the intensity 
along any radial line in the diffraction image of a circular aperture. 

In the diagram the curve of intensity along the major axis of the diffraction 
figure is called A, and that for the minor axis B. 

In the upper right-hand corner the curves are drawn together for the 



•03& 




5«IO- IX- 12- 13- 14- 15- 16- IT 



19- 20- 21- 2, 



■*>• 24- '•£• 26- M- 28' 29- &>■ 31- 52- 53' 54" 35- 36- 



Diagram 4. — Intensities of H, the Total Light Intensity along Axes of Symmetry of the 

Diffraction Figure. (Reduced to 1/5 linear of original.) 

A, Major axis and axis of symmetry of semicircle. 

B, Minor axis, perpendicular to axis of symmetry of semicircle. 

VOL. LXXXIII. — A, Z 
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central region of the image, the horizontal and vertical scales having as 
their respective units 8 = 1 and H = 0*05. For the outer regions of the 
image, it was found best to draw A and B separately and on an enlarged 
vertical scale, the new unit being H = 0*001, thus giving a magnification 
of 50, when compared with the central portion. 

The curve of B is drawn on the enlarged scale from 8 — 4 to 8 = 20, and 
A from 8 = 10 to 8 == 40. It will be seen that A has minima but no zeros. 

The curve B was not drawn beyond 8 = 20, as it is required only for 
comparison with A, and its last maximum shown is less than the last shown 
minimum of A. The values from which the curves have been drawn are 
given in Tables VIIIa and VIIIb for A and B respectively, and the methods 
by which the values were obtained are indicated. 



Table VIIIa. —Values of H (e = |tt). 



5. 



H. 



O'O 


1 -ooooo* 


0-4 


-98887* 


0-7 


0-96624* 


1-0 


-93212* 


1-4 


-87059* 


1-7 


-81422* 


2'0 


-75092* 


2*4 


-65887* 


2-7 


-58656* 


3*0 


-51358* 


3'3 


-44168* 


3'6 


-36920* 


3 -8318 


-32160* 


4-4 


-21091f 


5 -1356 


-10548f 


5-7 


-05182f 


6*0 


-03388f 


6-3 


-02119f 


6-7 


-01209f 


7 -0152 


-00986f 


7-5 


-01183f 


8-0 


-01825f 


8 -4172 


-02401f 


9"0 


-03056f 


9*5 


-03287f 



8. 



H. 



10-0 



10 
10 
10 



1734 

4 

9 

11-0 
11 -6199 
12-0 
12-5 
13*0 

13 -3238 

14-0 

14-1 

14 -7958 

15-0 

15-4 
15-9 

16-0 

16-470 

17-0 

18*0 

19-0 

19-5 



{ 



{ 

{ 
{ 



'03189$ 

(0 -0318)f 
-03087f 
-02912f 
-02828f 

(0'0231)f 
-01485f 

(0 -0107)f 
-00685f 

(0 -0047) f 
-00401f 
'00493$ 

(0'0049)f 
-00519f 
-00760f 
'00850$ 

(0 -0085)f 
-01004f 
-01041f 
-01056$ 

(0 -0104)f 
-01030f 
-00912$ 
'00560$ 
'00277$ 

(0 -00256) ] 



20 
21 

22 
22 
23 

24 

25 

26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 






5 























H. 







(0 


fo 

10 



















002801 
00359$ 
00490$ 
00504) t 



+ 



00485$ 
00348$ 
00348§ 
00195§ 
00137§ 
00187§ 
00268§ 
00296 § 
00238§ 
00149§ 
00097$ 
OOlll'l 
00161 § 
00192§ 
00172 § 
00119§ 
00070§ 
00074§ 
001045 



* Direct calculation (Taylor's theorem). 
$ Direct calculation (J's). 



( 



) Interpolated. 



t Mechanical integration. 
§ Semi- graphic method. 
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Table VIIIb. — Values of H (e = 0). 



5. 


H. 

1 


a. 


H. 


5. 


H. 


o-o 


1 -ooooo 


7-5 


'00130 


21-0 


-00027 


0-4 


-96066 


8-0 


-00344 


22*0 


-00011 


0-7 


-88356 


8 -4172 


-00416 


23'0 


o -ooooi 


1-0 


-77458 


9-0 


-00297 


24-0 


-00016 


1-4 


0*59910 


95 


-00115 


25-0 


o -oooio 


1*7 


-46203 


10 -1734 


o -ooooo 


26-0 


o -ooooo 


2'0 


-33261 


10-4 


-00011 


27-0 


o -oooio 


2*4 


0-18791 


10-9 


-00087 


28-0 


-00009 


2-7 


0-10700 


11 -6199 


-00161 


29-0 


-00002 


3-0 


-05109 


12-5 


-00070 


30-0 


-00006 


3-3 


-01789 


13 -3238 


o -ooooo 


31-0 


-00007 


3-6 


-00281 


14-1 


-00045 


32-0 


o -ooooo 


3 -8318 


-ooooo 


14 -7959 


-00078 


33-0 


-00004 


4 *4 


-00850 


15-4 


-00054 


34-0 


-00006 


5*1356 


-01750 


15-9 


-00019 


35-0 


o -ooooi 


5 7 


-01293 


16 -470 


o -ooooo 


36-0 


-00002 


6-0 


-00851 


17-0 


-00013 


37-0 


-00005 


6'3 


-00436 


18-0 


-00044 


38-0 


o -ooooi 


6-7 


-00081 


19*0 


-00012 


39-0 


o -ooooi 


7 '0152 


o -ooooo 


20-0 


-00005 

i 


40-0 


-00004 



Values obtained by direct calculation only. 



Method of Photographing the Diffraction Image. 

8. Although Scheiner and Hirayama have already published a photo- 
graph of the semicircular aperture diffraction figure, it seemed best to 
procure an independent picture of the fringes for comparison. Apparatus 
in use in the Observatory was employed as far as possible, and the only part 
specially made for the purpose was the aperture plate, of which a sketch is 
shown in Diagram 5. 

The base plate has a hole near the centre, into which, on the lower side, 
a chisel-shaped straight-edge is so fitted that the edge of the plane surface 
lies in a plane above the upper surface of the plate and parallel to it ; the 
edge is visible in the diagram in the centre of the hole. Over the upper 
surface of the plate, and a little above it, a carrier, with a circular hole in 
its centre and pivoted at one end, is moved parallel to the surface of the 
plate by a micrometer screw pushing the free end. The carrier is kept in 
good contact with the screw by a spiral spring. Into the circular hole of 
the carrier, a number of small circular discs may be fitted and fixed by two 
screws with washers, which are not shown. These discs have each a conical 
hole bored in the centre, the lower circular edge of which is in contact 
with the straight-edge already described. The two edges are kept in contact 
by a spring underneath the cross-bar shown near the micrometer screw. 
Finally the micrometer screw is fitted with a divided drum reading against 

z 2 
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Diagram 5. — Sketch of Apparatus employed. 
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a projecting scale. By this means the circular aperture is moved over the 
straight-edge until half of the circle is covered. The adjustment to the 
exact semicircle was made by means of a Hilger travelling microscope. This 
aperture was designed so that, after calibration of the micrometer screw, 
the diffraction fringes of any circular segment could be accurately found. 

The remainder of the apparatus used and the general arrangement is 
shown in Diagram 5. The equivalent focus of the telephoto-camera 
was about 19 feet. In consequence of the optical peculiarities of such a 
system it was necessary, if this focus was to be retained, to have the point 
source of light at a great distance, which, unfortunately, could not be done 
A collimator was accordingly introduced ; and, as the image of the pinhole 
source of light was magnified in the ratio of the foci of camera and 
collimator, a small image of the pinhole produced by an eye-piece was used, 
the only disadvantage being a largely increased exposure. The exposure 
required to obtain the photograph reproduced (Plate 3) was three hours. 
The greatest difficulty met with arose from the object-glass of the telescope 
having a number of waves of unequal refraction in the glass, and the parts so 
affected had to be most carefully avoided, as otherwise the diffraction image 
would have been adversely affected. 
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Comparison of Results, 

9. We have now a contour diagram, two photographs and a diagram of 
the intensity along the principal axes of the diffraction image. How do 
these different results agree together ? 

We have first to compare the results given for maxima and minima by the 
photographs and intensity curves. It was not possible to measure accurately 
the optical constants of the photographic apparatus used, so that the only 
available method of determining the reduction factor of my own photograph 
was to measure the dimensions along the minor axis— where the diffraction 
image is the same as that of a circular aperture. We have the following 
results : — 

Table IX. — Major Axis. 





By intensity 
diagram. 


By enlargement 

of own 

photograph. 


By enlargement of 

Scheiner and Hirayama's 

plate. 


2nd min 

2nd max 

3rd min. , 


14-2 
19*4 
26*8 
32'0 
38*9 
44-9 


13-9 
19-4 

31-8 

44-2 


14-3 
19-9 
26-8 
32*6 

45-1 



These figures give the diameter of the image along the major axis bounded 
by the particular points of the image mentioned, and measured in terms of 
the adopted unit S. The position and intensities of the maxima and minima 
along both axes are given below : — 



Table X. 



Major axis. 



5. 



Centre 

1st min 

1st max 

2nd min 

2nd max 

3rd min , 

3rd max 

4th min 

4th max 

5th min 

5th max t . 

6th min 



O'O 
7-1 
9-7 
13-4 
16*0 
19-45 
22-45 
25-95 
28-85 
32 -25 
35-15 
38*50 



H. 



Minor axis. 



5. 



H. 



1-000 


o-o 


0-009 


3*832 


0-033 


5-136 


-00397 


7 -015 


-01057 


8-417 


-00257 


10 -173 


-00503 


11 -620 


-00137 


13 -324 


-00296 


14 -796 


-00095 


16 -470 


-00193 


17 -960 


-00068 


19 -616 



1 -ooooo 
o -ooooo 

-01750 

o -ooooo 

-00416 

-ooooo 

-00161 

-ooooo 

-00078 

o -ooooo 

-00043 

-ooooo 
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The values for the major axis were obtained from the intensity curve 
shown in Diagram 4, while those for the minor axis are, of course, the 
values for the diffraction image due to a circle. There is an interesting point 
about these values. In the case of the minor axis (or circle) each maximum 
lies nearer to the minimum of its own order than to that of the next higher 
order, but this tendency of the maxima decreases as the order increases, and 
ultimately vanishes. In the case of the major axis, the same tendency 
of the maxima is present ; but it appears that the tendency, although 
becoming smaller, does so far less rapidly. 

Turning now to the contour diagram (fig. 2), and comparing it with the 
photograph (Plate 3), the agreement is perfect ; and while the photograph 
makes reference to the contour diagram easier, the contour diagram emphasises 
and extends the details of the photograph. Erom the contour diagram we see 
that the central spot is not really an ellipse, but only an approximation to 
one. Calling the five contour lines shown in the central spot A, B, C, D, E, 
counting from the centre outwards, we find, on measuring their diameters 
along the axes of the figure : — 

Table XI. 





Major axis. 


Minor axis. 


Ratio. 


A 


2 '436 


1-296 


1-88 


B 


6 -114 


3 -234 


1-90 





8*814 


4 -720 


. 1-87 


D 


11 -446 


6-020 


1-90 


E 


12 '686 


6*550 


1-94 



We see that there is, on the whole, an increase in the ■ ratio of the 
diameters, so that the shape of the central spot varies with the exposure 
when taking a photograph, and also with the intensity of the light when 
observed visually. The deviations from the steady increase in the ratio 
which are shown by B and C may be quite reasonably explained by the 
influence of the first bright ray, which should increase the ratio very slightly 
in the case of B, and of the first dark ray, which appears to decrease the 
ratio somewhat in the case of C. In connection with this, the way the 
contour line we have just called E is pulled out of shape is very significant. 

10. General description of the diffraction image. 

The image is, as has been already mentioned, symmetrical about two- 
perpendicular axes. There is a central approximately oval spot, which 
is roughly a little less than twice as long as it is broad. This central spot 
is bisected along its major axis by a long bright ray, which fades, away as it 
recedes from the centre. Around the spot and following its outline are a 
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series of alternate dark and bright rings, which, at the minor axis, are 
identical with those of the diffraction image of a circular aperture, and, 
along the major axis, terminate at an angle in a series of bright spots. 
Further, these rings are crossed and modified by a series of alternate bright 
and dark rays of approximately hyperbolic shape, having the major axis of 
the figure for their geometric axis, which they cut at its maxima and 
minima, with the exception that the first one of the series passes through a 
point in the central bright spot. The whole diffraction image presents an 
exceedingly beautiful appearance, which is, unfortunately, only very 
inadequately represented by photography, but which is very easily studied 
with the double-star apparatus when monochromatic light is used and the 
star is examined in the equatorial through the heliometer lens. 

Conclusions. 

(i) The complex diffraction fringes due to a semicircular aperture have 
been for the first time theoretically determined, and actual measurements on 
the fringes themselves shown to be in close accordance with the calculated 
results. The absolute agreement of the computed contours with the 
photographs forms another item in the evidence, showing the old 
undulatory theory of light suffices to describe very accurately diffraction 
phenomena. 

(ii) A problem which has attracted the attention of more than one 
astronomer and mathematician but had not yet been overcome, has been 
solved, and the exact nature of heliometer images numerically determined. 

(iii) The solution of the problem has illustrated the limits to the practical 
utility of Bessel's functions ; the mechanical integration of 24 indefinite and 
17 definite integrals has indicated the manner in which graphical methods 
can be used to supplement numerical analysis, and has contributed to the 
solution of a physical and astronomical problem of considerable interest, 
which had hitherto proved unattainable by the usual methods of evaluation. 
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